Yang-Lee and Fisher zeros generalized on some far-from-equilibrium systems. 
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A generalization of the Yang-Lee and Fisher zeros on far-from-equilibrium systems coupled with 
. . . two thermal baths is proposed. The Yang-Lee zeros were obtained for minimal models which ex- 

hibit complicated behavior in the context of the partition function zeros and provide an analitycal 

■ treatment. This type of models may be considered as a simplest one and analogous to Ising model 
^ ' for equilibrium. The obtained distributions of generalized Yang-Lee zeros show nontrivial behavior 

for these simple models. 
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, More than five decades ago Yang and Lee proposed an approach to clarify how singularities of the thermodynamic 
' functions appears within the canonical ensemble [J . They considered liquid-gas transition and wrote down partition 
function in grand-canonical description as a polynomial with respect to fugacity. The main idea was an analysis of the 
fugacity in complex plane. Complex roots of this polynomial with respect to fugacity cross the real positive semi-axis 
O in the singularity points in the thermodynamic limit and single out the points of the phase transition. Further, these 
^ roots in the complex fugacity plane are called Yang-Lee zeros. After Yang and Lee pioneer work many papers are 
S ' devoted to this description of the phase transitions. 

^Jj , The approach of Yang-Lee zeros was used to describe spin system where complex fugacity in grand-canonical 

■ ensemble was replaced by the complex value of e~^^ with h as a magnetic field. It was shown that the density of the 
c/2 . Yang-Lee zeros is in a close relation with the critical exponents describing the phase transition and may be used as a 

measure of the strength of the transition. Besides the Yang-Lee zeros so-called Fisher zeros were invented in complex 
temperature plane [4]. Their properties are more dependent on the particular system. For the historical details and 
y ■ useful references on Yang-Lee and Fisher zeros see review Here we describe only some general aspects. 
K ' It was shown by Binek et a/.fi*] that the Yang-Lee zeros are not fully and only a theoretical concept. They proposed 
(— ( . an experimental way to measure the Yang-Lee edge singularity exponents from the isothermal magnetization data in 
Q ' 2D Ising ferromagnet. 

, In summary, the partition function zeros are rather well investigated in case of the systems in equilibrium and 
provides us a clue in understanding of the phase transitions. One can measure some quantities which are directly 
I [ related to the partition function zeros and get some experimental evidence. The question is how generalize an approach 
^ ■ of Yang-Lee zeros to nonequilibrium systems and when such a possibility exists. 

^\ A standart approach to describe the behavior of the nonequilibrium system is to derive master equation and 

C ■ determine the rates of the processes in it. After a long time, the system may settle in a nonequilibrium steady- 
state. The exact meaning of the "long time" depends on the rates of processes in it. According to initial state 
and values of the parameters describing the system there may be different steady states and the transition from one 

■ of them into one another may be considered as an analogy to equilibrium phase transition. The one possibility to 
J generalize the partition function approach is a consideration of the nonequilibrium models which allow the transfer 

25 ' matrix description. As it is proposed by Arndt 0] some nonequilibrium models which are useful for the description 
of the driven diffusive systems, traffic flow, biological transport and other processes give us a possibility to invent the 
partition function and its zeros. In these models we can express the stationary probability distribution as a trace over 
some algebra elements, formulate the concept of partition zeros and show the direct relation between steady-state 
^ I , transition and the distribution of the partition function zeros. However, it is possible to derive the distribution of 
the Yang-Lee zeros without having the exact expression for the partition function. One can consider the zeros of a 
steady-state normalization factor in the complex plane of the transition rates @ . 

In this work one another possibility to generalize the the partition function zeros is proposed. Let us consider 
thermodynamics of two temperature systems with different time-scales and Hamiltonian H{a, s) depending on fast 
variables a coupled with thermal bath at temperature and slow variables s coupled with bath at temperature T^. 
The system may be considered as a minimal model having easily controlled nonequilibrium properties . The one of 
the real examples may be NMR/ESR physics, where one of the baths is realized by weak dipole interactions and the 
second one is a lattice temperature 8] . 

In such a system one can expect that after some long time system settles in a nonequlibrium steady-state with 
heat currents between the baths. This state is called steady adiabatic 0. Although the stationary distribution is 
far from Gibbsian, one can derive the Gibbs like corresponding stationary distribution. As it is shown in [lo[ this 
type of nonequilibrium systems has complicated and unusual behavior from equilibrium point of view. There may be 
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nonequilibrium phase transitions with no equihbrium counterpart. Also, it is possible for latent heat to be negative in 
presence of conflicting interactions [l0| . We expect that the properly generalized partition function zeros may violate 
some usual equilibrium statements too. 

According to 10] and taking into account a huge difference of rates between fast and slow variables one can write 
down the conditional probability 

^('^l^) = -7^,^~^"^"''\Z{s) = rr,e-^^(^^^). (1) 
Z{s) 

Here we assume that on the times relevant to cr, s does not change and the conditional probability Eq. ([1]) is Gibbsian. 
The next step is to derive P{s). This may be done by averaging the force dsH{a, s) acting on s over the distribution 
P{a\s) 

P{s)^^Z^{s),Z = Trse-^F{s), (2) 

where F{s) = —T^lnZ{s) is the conditional free energy and the common probability is P(cr, s) = P{s)P{(7\s). 
It is possible to obtain the steady state distribution by minimizing the free energy 

F = -TJnZ. (3) 

Eq. ^ shows the way to define the partition function zeros. It may be done like the equilibrium one, after replacing 
the equilibrium partition function by Z in common equilibrium description. As it is in equilibrium, one may suppose 
that there is a possibility of the far-from-equilibrium phase transition between steady states in thermodynamic limit 
[Tot , when the probability distribution is quasi- Gibbsian ^ (as it is for Gibbsian). The example of phase transition 
is demonstrated in [lol |. However, this simple example is still mathematically rather complicated. The motivation of 
usefulness of the partition function zeros and their application for the indication of the phase transition points in case 
of the quasi-Gibbsian distribution is the same as it is in equilibrium. The main purpose of this work is not to solve 
complicated models demonstrating phase transitions and rather difficult for exact solution. Our aim is to investigate 
the properties of the partition function zeros for the models providing analytical solution and compare the results 
with equilibrium analogs. 

Let us discuss the simplest models which may be treated analytically and define the partition function zeros on those 
examples. As such simple systems one can take the one-dimensional Ising model with nearest-neighbor interactions 
and couplings or magnetic fields as fast variables 

N N 

H ^ -^JiSiSi+i -^hiSi. (4) 

i=l i=l 

For the simplicity, let the couplings Ji are the fast variables and take values ±1. In that case one may write down for 
Z{s) the following expression 

Z{s) = Tr^.n}e-^" = 2^\{coM'^)e^'\ (5) 

i 

where Tj is a temperature of the thermal bath coupled to fast variables. From Eq. ([5]) we get for Z 

Z = Tr[s,}Z{s)Ti ^ TrV^, (6) 

where the transfer-matrix 

y=[2cosh(£|±i)]5fe*^, (7) 

is introduced. Here we assume the periodic boundary conditions and that the system consists of N spins. The 
transfer-matrix V is symmetric and one can write Z with respect to eigenvalues of V, as it is done in equilibrium case 
for one-dimensional models 



The condition Z = leads to 



Z = \^ + X^. (8) 
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One of the examples of this procedure and its more detailed description in equilibrium case may be found in [llj . 

As a generalization of the equilibrium Yang-Lee zeros one can consider solution of Eq. ([9|) with respect to /x = e . 
However, the model with fast random couplings is too trivial. After taking into account symmetry of the function 
cosh and some simple algebra one gets the set of the Yang-Lee zeros as /i = ±i. A similar set of values fi = (±1)-^'*/-^= 

h 

would be obtained if we generalize Yang-Lee zeros as /x = e . 

In case of fast magnetic fields hi = zth the corresponding transfer-matrix is 



V 



2 cosh 



h{si + Si+i) 



2Th 

It is easy to obtain for the eigenvalues of Eq. (fTO)) the following expression 



(10) 



Al,2 — 



2 cosh — 



h 



eT7±2~e~*^. (11) 



Here the natural and most simplest choice of the generalized Yang-Lee zeros is solution of Eq. © with respect to 

h h h 

e "^h . For both e and e the general assumption is a complex magnetic field and the corresponding temperature 
is chosen according to details of model. It is easy to obtain the second type of the Yang-Lee zeros if the first one is 
known and vice versa. 

h T\ _J_ 

Denoting by /i = e and a = 2 e we obtain 

+ ^e^ + a = {[{^l + ^ ^ ay^. (12) 

Using condition on ([5]) one can derive a finite set of independent equations corresponding to all possible values of 
if. These equations for generalized Yang-Lee zeros are transcendental in general. One can consider a simple choice 
of rational ratio ^ or for more simplicity we assume it to be a natural number. It is possible for the natural ratio 

^ to solve the set of equations numerically for even a large numbers oi N 10^. The obtained distribution may 
be considered as to be close to the thermodynamic limit. However, the fast magnetic fields model has an obvious 
analytical solution in thermodynamic limit. Near the thermodynamic limit one can treat ip as a, continuous parameter, 
running over the interval (0, 27r) due to relation N oo. One can derive from Eq. the following relation 

[(. + ,.-.)]* ^2*e-*f:^. (13) 

As an example, if the ratio is = 2 one gets two simple quadratic equations and the corresponding solutions are 



Ml,2,3,4 = ^ , (14) 

where c — 2e~^ . The corresponding distribution is represented in Fig. [T] Four solutions of Eq. ([T3]) form a 

complicated picture. This behavior is rather different from equilibrium one, when for a simple Ising-type model one 
gets partition function zeros lying on a unit circle. Also, it is easy to detect from Eq. that for ip ^ the both 
four solutions are /ii ■2..;^.4 — > oo. In equilibrium case the usual condition on Yang-Lee edge singularity points is f = 
(see for examples [ij, [iJl ). The second step one can do is to consider model with fast variables ji = ±7, slow 
variables s; = ±1 and Hamiltonian as 

N N 

H = -^jiS^Si+i -^jiSi. (15) 
1=1 1=1 

The corresponding transfer-matrix is 

V=[2cosh(2f^-,^ifii^)]??. (16) 

It is straightforward to define the generalized Yang-Lee zeros for this model, as partition function zeros with respect 
to n = e'*^'^^ . After writing down characteristic equation for transfer-matrix and obtaining the eigenvalues it is simple 
to get the Yang-Lee zeros due to Eq. After some numerical calculation within Mathematica we obtained a highly 
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FIG. 1: The distribution of generalized Yang-Lee zeros in complex fj, — e'^h plane, where J — 0.5 and Th = 10, Ts — 5. 
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FIG. 2: The distribution of generalized Yang-Lee zeros in complex fi = e'^i plane, where N — 10^ and T^/Th = 2. 



nontrivial distribution for such a simple model (see Fig. [5] ) . The behavior of the zeros (Fig. ^ also violates the circle 
theorem. 

In order to generalize the Fisher zeros one can introduce two types of the Fisher zeros according to two temperatures. 
In equilibrium limit, when the ratio ^ is equal to one those two types of the Fisher zeros coincide. Simple algebra 
shows that the same pointless set we obtain for the Fisher zeros in case of fast couplings. For the model with fast 
magnetic fields one need to solve more complicated equations compared with those for Yang-Lee zeros. 

Summarizing, the basic approach to generalize the partition function zeros for non-equilibrium systems coupled to 
two thermal baths is proposed. The described zeros show nontrivial behavior and need to more complete investigation. 
This approach of the generalized partition function zeros may be used to analyze phase transitions by numerical and 
analytical methods as it is done in equilibrium case. As author supposes, it is interesting to discuss the properties of 
nonequilibrium Yang-Lee and Fisher zeros in more details and for more complicated systems. 

Author is grateful to A. E. AUahverdyan for bringing attention to this field and N. Ananikian for useful discussions. 
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